Algebraic study of drifting spiral waves.

This two-dimensional study is motivated by cardiac electrophysiology, and focuses on rotating spiral waves in reaction-diffusion (RD) models. Here we deal with a spiral's translational drift under a constant externally imposed gradient G. A long-standing problem may be stated as follows: Given the dimensionless drift velocity V/G, find its nontrivial direction angle Γ relative to G. A deductive algebraic treatment yields a solution, cosΓ=-V/G. Three features are worth noting: the combination of algebraic and RD contexts; a somewhat extensive derivation contrasting with a compact result; and the generality due to the absence of reaction details in the formula. Agreement with a computational database is good to fair, if spirals of very low density are excluded.

Generalized minimal principle for rotor filaments.

To a reaction-diffusion medium with an inhomogeneous anisotropic diffusion tensor D, we add a fourth spatial dimension such that the determinant of the diffusion tensor is constant in four dimensions. We propose a generalized minimal principle for rotor filaments, stating that the scroll wave filament strives to minimize its surface area in the higher-dimensional space. As a consequence, stationary scroll wave filaments in the original 3D medium are geodesic curves with respect to the metric tensor G=det(D)D(-1). The theory is confirmed by numerical simulations for positive and negative filament tension and a model with a non-stationary spiral core. We conclude that filaments in cardiac tissue with positive tension preferentially reside or anchor in regions where cardiac cells are less interconnected, such as portions of the cardiac wall with a large number of cleavage planes.

Rotation of a scroll wave reveals the metric tensor for the associated geodesic filament.

This study considers a reaction-diffusion medium that supports a rotating scroll wave with stationary filament and diffusivity matrix D. The filament is a geodesic curve in three-space. The metric tensor g associated with that geodesic is obtained in terms of D, resulting in g=(detD)D-1. Here detD (the determinant of D) is not necessarily constant in space. The detD factor, although strongly indicated previously, has not yet been derived deductively, as it is here. We take the opportunity to outline, starting from the reaction-diffusion equations, a sequence of mathematical steps leading to both geodesic and metric.

Refraction of scroll-wave filaments at the boundary between two reaction-diffusion media.

We explore the shape and the dynamics of scroll-wave filaments in excitable media with an abruptly changing diffusion tensor, important for cardiac applications. We show that, similar to a beam of light, the filament refracts at the boundary separating domains with different diffusion. We derive the laws of filament refraction and test their validity in computational experiments. We discovered that at small angles to the interface, the filament can become unstable and develop oscillations. The nature of the observed instabilities, as well as overall theoretical and experimental significance of the findings, is discussed.

Probing field-induced tissue polarization using transillumination fluorescent imaging.

Despite major successes of biophysical theories in predicting the effects of electrical shocks within the heart, recent optical mapping studies have revealed two major discrepancies between theory and experiment: 1), the presence of negative bulk polarization recorded during strong shocks; and 2), the unexpectedly small surface polarization under shock electrodes. There is little consensus as to whether these differences result from deficiencies of experimental techniques, artifacts of tissue damage, or deficiencies of existing theories. Here, we take advantage of recently developed near-infrared voltage-sensitive dyes and transillumination optical imaging to perform, for the first time that we know of, noninvasive probing of field effects deep inside the intact ventricular wall. This technique removes some of the limitations encountered in previous experimental studies. We explicitly demonstrate that deep inside intact myocardial tissue preparations, strong electrical shocks do produce considerable negative bulk polarization previously inferred from surface recordings. We also demonstrate that near-threshold diastolic field stimulation produces activation of deep myocardial layers 2-6 mm away from the cathodal surface, contrary to theory. Using bidomain simulations we explore factors that may improve the agreement between theory and experiment. We show that the inclusion of negative asymmetric current can qualitatively explain negative bulk polarization in a discontinuous bidomain model.

Frustrated drift of an anchored scroll-wave filament and the geodesic principle.

We investigate anchored scroll-wave filaments in an excitable medium whose diffusivity matrix, including its determinant, is spatially nonuniform. The study is motivated by cardiological applications where scroll-wave behavior in the presence of diffusivity gradients is believed to play an important role in the development of severe arrhythmias. A diffusivity gradient is expected to make the filament drift, unless drift is prevented ("frustrated") by anchoring to localized defects in the propagation medium. The resulting stationary filament is a geodesic curve, as demonstrated here in the case of a nonzero but constant gradient. That is, the diffusivity matrix has a determinant that varies in space, in contrast to what was assumed in earlier work. Here, we show that the filament shape results from a metric tensor of the form (det D)D{-1} , where D is the diffusivity tensor. The filament's shape is solely determined by the diffusivity tensor and is independent of the equation's reaction terms. We derive the analytic solution for the filament and determine conditions for the existence of that solution. The theory is in excellent agreement with numerical simulations.

Monitoring intramyocardial reentry using alternating transillumination.

Intramyocardial reentry is implicated as a primary cause of the most deadly cardiac arrhythmias known as polymorphic ventricular tachycardia and ventricular fibrillation. However, the mechanisms involved in the triggering of such reentry and controlling its subsequent dynamics remain poorly understood. One of the major obstacles has been a lack of adequate tools that would enable 3D imaging of electrical excitation and reentry inside thick ventricular wall. Here, we present a new experimental technique, termed alternating transillumination (AT), aimed at filling this gap. The AT technique utilizes a recently synthesized near-infrared fluorescent voltage-sensitive dye, DI-4-ANBDQBS. We apply AT to study the dynamics of reentry during shock-induced polymorphic ventricular tachycardia in pig myocardium.

Multiplicative optical tomography of cardiac electrical activity.

Cardiac electrical activity can be mapped today through the response of voltage-sensitive dyes; but poor transparency of muscle tissue has enforced shallow-depth imaging. We present a three-dimensional (3D) reconstruction method for electrical activity deep inside the myocardial wall. Our approach is nonlinear and differs substantially from standard diffusive optical tomography. It does not require matrix inversion, data regularization or a priori information concerning the original object. Opposite sides of a slab-shaped preparation are scanned in parallel by detection and illumination points with a constant vector offset between illumination and detection axes (biaxial scanning). Scanning is performed in two perpendicular directions. In each direction, a pair of 2D images is obtained under offsets of opposite signs. These two pairs are the input for a multiplicative reconstruction algorithm, whose output is a 3D image. The overall procedure was successfully tested on computer-generated sources that include points, lines and hemispheres, patterned after actual electrophysiological excitations. The algorithm is computationally efficient and stable with respect to varying noise levels in the raw data.

Asymmetric bound states of spiral pairs in excitable media.

We present a stable regime of asymmetric bound states for spiral pairs in a generic numerical model of a homogeneous excitable medium. In this regime, one spiral tip (slave) rotates around the other (master). Master-slave dynamics occur for both same-chirality and opposite-chirality spiral pairs in a range of parameters and initial conditions. We study the dependency of master-slave characteristics on the medium's excitation threshold and present a phenomenological model that accounts for the qualitative properties of master-slave dynamics.

Synthesis of voltage-sensitive fluorescence signals from three-dimensional myocardial activation patterns.

Voltage-sensitive fluorescent dyes are commonly used to measure cardiac electrical activity. Recent studies indicate, however, that optical action potentials (OAPs) recorded from the myocardial surface originate from a widely distributed volume beneath the surface and may contain useful information regarding intramural activation. The first step toward obtaining this information is to predict OAPs from known patterns of three-dimensional (3-D) electrical activity. To achieve this goal, we developed a two-stage model in which the output of a 3-D ionic model of electrical excitation serves as the input to an optical model of light scattering and absorption inside heart tissue. The two-stage model permits unique optical signatures to be obtained for given 3-D patterns of electrical activity for direct comparison with experimental data, thus yielding information about intramural electrical activity. To illustrate applications of the model, we simulated surface fluorescence signals produced by 3-D electrical activity during epicardial and endocardial pacing. We discovered that OAP upstroke morphology was highly sensitive to the transmural component of wave front velocity and could be used to predict wave front orientation with respect to the surface. These findings demonstrate the potential of the model for obtaining useful 3-D information about intramural propagation.

Minimal principle for rotor filaments.

Three-dimensional rotors, or scroll waves, provide essential insight into the activity of excitable media. They also are a suspected cause in the formation and maintenance of ventricular fibrillation, whose lethality is well known. It is therefore of considerable interest to find out what configurations can be adopted by such pathologies. A scroll's behavior is embodied in its organizing center or filament, a largely quiescent tube about which the scroll rotates. Predicting filament shape has normally required computer-intensive simulations of the whole scroll in time. We have found a fast and robust principle that yields the prediction for stationary filaments on a purely geometrical basis, blind to the reaction parameters of the medium. The procedure is to calculate the filament shape as a minimal path. We work in singly diffusive media whose diffusivity tensor--and no other feature--varies spatially. Mathematical and numerical evidence is presented for the proposition that a stable filament is a geodesic in a three-dimensional space whose metric is given by the inverse diffusivity tensor of the medium. Away from the boundaries, a stable filament is unaffected by the reaction parameters. The algorithmic aspects of this work are subsidiary to our main purpose of drawing attention to the universal and unexpectedly exact fit of an elementary geodesic principle within reaction-diffusion theories.